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Algorithmic Cooling (AC) of spins applies entropy manipulation algorithms in open spin-systems 
in order to cool spins far beyond Shannon's entropy bound. AC of nuclear spins was demonstrated 
experimentally, and may contribute to nuclear magnetic resonance (NMR) spectroscopy. Several 
cooling algorithms were suggested in recent years, including practicable algorithmic cooling (PAC) 
and exhaustive AC. Practicable algorithms have simple implementations, yet their level of cooling 
is far from optimal; Exhaustive algorithms, on the other hand, cool much better, and some even 
reach (asymptotically) an optimal level of cooling, but they are not practicable. 

We introduce here semi-optimal practicable AC (SOPAC), wherein few cycles (typically 2-6) are 
performed at each recursive level. Two classes of SOPAC algorithms are proposed and analyzed. 
Both attain cooling levels significantly better than PAC, and are much more efficient than the 
exhaustive algorithms. The new algorithms are shown to bridge the gap between PAC and exhaustive 
AC. In addition, we calculated the number of spins required by SOPAC in order to purify qubits for 
quantum computation. As few as 12 and 7 spins are required (in an ideal scenario) to yield a mildly 
pure spin (60% polarized) from initial polarizations of 1% and 10%, respectively. In the latter case, 
about five more spins are sufficient to produce a highly pure spin (99.99% polarized), which could 
be relevant for fault-tolerant quantum computing. 



I. INTRODUCTION 

We consider here spin-half nuclei, quantum bits 
(qubits), whose computation-basis states align in a mag- 
netic field either in the direction of the field (|f) = |0)), 
or in the opposite direction (||) = |1)). Several such bits 
(of a single molecule) represent a binary string, or a reg- 
ister. A macroscopic number of such registers/molecules 
can be manipulated in parallel, as done, for instance, in 
nuclear magnetic resonance (NMR). From the perspec- 
tive of NMR quantum computation (NMRQC) [1] (for a 
recent review see [2] ) , the spectrometer that monitors and 
manipulates these qubits/spins can be considered a sim- 
ple "quantum computing" device (see for instance Refs. 
[3] and [4] and references therein). The operations (gates, 
measurements) are performed in parallel on many regis- 
ters. 

The probabilities of a spin to be up or down are 
p^ = i^p 1 and p^ = where £o is called the po- 

larization bias [5]. At equilibrium at room temperature, 
nuclear-spin polarization biases are very small; In com- 
mon NMR spectrometers £o is at most ~ 1CP 5 . While 
the polarization bias may be increased by physical cool- 
ing of the environment, this approach is very limited for 
liquid-state NMR, especially for in- vivo spectroscopy. To 
improve polarization, including selective enhancement, 
various "effective cooling" methods increase the polar- 
ization of one or more spins in a transient manner; The 
cooled spins re-heat to their equilibrium polarization as a 
result of thermalization with the environment, a process 
which has a characteristic time of Ti. 

An interesting effective cooling method is the re- 
versible (or unitary) entropy manipulation of S0rensen 
and Schulman-Vazirani [6, 7], in which data compression 



tools and algorithms are used to compress entropy onto 
some spins in order to cool others. Such methods are 
limited by Shannon's entropy bound, stating that the 
total entropy of a closed system cannot decrease. The 
reversible cooling of Schulman-Vazirani was suggested as 
a new method to present scalable NMRQC devices. The 
same is true for the algorithmic cooling approach [8] that 
is considered here. 

Algorithmic Cooling (AC) combines the reversible 
effective-cooling described above with thermalization of 
spins at different rates (Ti values) [8]. AC employs slow- 
thermalizing spins (computation spins) and rapidly ther- 
malizing spins (reset spins). Alternation of data com- 
pression steps that put high entropy (heat) on reset spins, 
with thermalization of hot reset spins (to remove their ex- 
cess entropy) can reduce the total entropy of suitable spin 
systems far beyond Shannon's bound. Let us describe in 
detail the three basic operations of AC. 

1. RPC. Reversible Polarization Compression 

steps redistribute the entropy in the system so that 
some computation spins are cooled while other com- 
putation spins become hotter. 

2. SWAP / POLARIZATION TRANSFER. 

Controlled reversible interactions allow the hotter 
computation spins to adiabatically lose their en- 
tropy to a set of designated reset spins. Alterna- 
tively, we say that the higher polarization of the 
reset spins is transferred onto these hotter compu- 
tation spins. 

3. WAIT / RESET. The reset spins rapidly ther- 
malize (much faster than the computation spins), 
conveying their excess entropy to the environment, 



while the computation spins remain colder, so that 
the entire system is cooled. 

In some algorithms, the reset spins can be used directly 
in the compression step [9, 10]. These algorithms have 
already led to some experimental suggestions and imple- 
mentations [11-16]. 

An efficient and experimentally feasible AC technique, 
termed "practicable algorithmic cooling (PAC)" [9], com- 
bines polarization transfer (PT), reset, and RPC among 
three spins called 3-bit- compression (SB- Comp). The 
subroutine 3B-Comp may be implemented in several 
ways. One implementation is as follows [9]: 
Exchange the states 1 100) «-» 1 01 1) . 

Leave the other six states (\000), |001), etc.) unchanged. 
If 3B-Comp is applied to three spins {C,B,A} with bi- 
ases £o <C 1, then spin C will acquire the bias [9]: 
e' = ^j 1 . PAC1 and PAC2 consider purification levels, 
where spins are cooled by a factor of 1.5 at each suc- 
cessive level; At the (j + l) st purification level, ej+i = 
1.5sj = 1.5 j+1 e [17]. ' 

The practicable algorithm PAC2 [9] uses any odd num- 
ber of spins, n = 2 J + 1: one reset spin and 2 J compu- 
tation spins (see Appendix A for a formal description 
of PAC2). PAC2 cools the spins such that the coldest 
spin can, ideally [18], attain a bias of (3/2) J . When 
all spins are cooled, the following biases are reached: 
{. . . , 5^, 3f , 3f , 2± 2±, if, 1±, 1, 1} . This proves an ex- 
ponential advantage over the best possible reversible 
cooling techniques (e.g., of Refs. [6] and [7]), which are 
limited to a cooling factor of y/n. As PAC2 can be ap- 
plied to as few as 3 spins ( J = 1) to 9 spins ( J = 4), using 
NMRQC tools, it is potentially suitable for near future 
applications [12]. PAC is simple, as all compressions are 
applied to three spins and use the same subroutine, 3B- 
Comp. 

PAC2 always applies compression steps (3B-Comp) to 
three identical biases. In general, for three spins with 
biases, (ec,eb,ea), where all biases <C 1, 3-bit compres- 
sion would attribute spin C with the bias [15, 19, 20]: 

i _ e c + £b + £a m 
£c ~ 2 ' ( j 

Using Eq. 1, several algorithms were designed, including 
PAC3, which always applies 3B-Comp to bias configura- 
tions of the form {sj, Ej, £j-i} (see formal description in 
Appendix A). When all spins are cooled, the following 
biases are obtained: {. . . , 9^, 7, 5g, 3|, 2|, 2, 1^, 1, 1} . 

The framework of PAC was extended to include mul- 
tiple cycles at each recursive level. Consider the simplest 
case, whereby a three-spin system (C,B,A), with equal 
bias £o ; is cooled by repeating the following procedure m 
times: 

1. 3B-Comp on (C, B, A), increasing the bias of C. 

2. reset the biases of A and B back to £ . 

The bias of spin C increases to 2eo(l — 2~ m ). Thus 
the biases {2,1,1} are asymptotically obtained for the 



three spins, as noted first by Fernandez [19]. This vari- 
ant, where m 3> 1, is the first exhaustive spin-cooling 
algorithm. By generalizing this "Fernandez algorithm" 
to more spins, we obtained [15, 21, 22] a bias configu- 
ration that asymptotically approaches the Fibonacci se- 
ries; when all spins are cooled, the following biases are 
reached: {. . . , 34, 21, 13, 8, 5, 3, 2, 1, 1} . Algorithms based 
on the 3B-Comp were recently reviewed [23]. 

Following the Fibonacci algorithm, a related algorithm 
was described, which involves compression on four ad- 
jacent spins [15, 22]; The tribonacci algorithm reaches 
biases in accord with the tribonacci series (also called 
3-step Fibonacci series) {...,81,44,24,13,7,4,2,1,1}, 
where each term is the sum of the three previous terms. 
Similarly, one can obtain general fc-step Fibonacci se- 
ries, where each term is the sum of the previous k el- 
ements, respectively. Summing over all previous terms 
produces the geometric sequence with powers of two: 
{...,128,64,32,16,8,4,2,1,1}. The AC that reaches 
this series was termed [15, 22], "all-bonacci AC". We 
conclude that, with the same number of spins, exhaus- 
tive AC reaches a greater degree of cooling than PAC 
algorithms. For example, the cooling factor achieved by 
Fibonacci with 9 spins is T§ = 34. The corresponding 
enhancement for PAC 2 is (3/2) 4 w 5, and for PAC3 it is 
about 9^. 

The Partner Pairing Algorithm (PPA) was shown to 
achieve a superior bias than all previous cooling algo- 
rithms, and was proven to be optimal [10, 21]. The all- 
bonacci algorithm appears to produce the same cooling 
as the PPA: {...,16,8,4,2,1,1} (this was verified nu- 
merically for small biases [15, 22]). In all-bonacci, with 
n spins, the coldest spin approaches a bias of £/ = 2 n ~ 2 £ 
(as long as £/ <C 1). For the PPA, a calculation of the 
cooling that yields 2™~ 2 £o was not performed, yet it was 
proven that '2"- 2 £ < s f < 2™- 1 £ . See [10, 21] for 
a proof and see [15, 22] for calculations in the case of 
Ef « 1. 

The focus of this paper is on simple algorithms that 
are semi-optimal and also practicable, owing to their ex- 
clusive reliance on simple logic gates. Such algorithms 
could be implemented in the lab if proper physical con- 
ditions (e.g., Ti ratios) are achieved, and thus they 
could have practical implications in the near future. In 
contrast, the optimal algorithms mentioned above, the 
PPA and the all-bonacci, do not belong to this category. 
The all-bonacci requires an unreasonable number of reset 
steps [24], and the permutations required in the sorting 
step of the PPA were not yet translated into simple 2-spin 
and 3-spin gates [25]. 

The next section describes rnPAC, a new semi-optimal 
algorithm based on PAC, and its special variant, mPAC, 
which cool n = 2 J + 1 spins such that the coldest spin 
is asymptotically cooled by a factor of (2 — 2~ m ) J . The 
optimal algorithms, the PPA and all-bonacci, need about 
half the spins (more precisely, J+2 spins) in order to cool 
(asymptotically) the coldest spin to 2 J ; we thus see that 
mPAC is semi-optimal in the sense that it needs twice 
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as much spins to reach nearly the same optimal bias. 
Section III discusses simple variants of the Fibonacci al- 
gorithm: (^-Fibonacci [21, 23], and mFib, which fixes the 
number of cycles at each recursive level. Section IV com- 
pares the new algorithms, mPAC and mFib, to previous 
cooling algorithms, including the PPA, Fibonacci, all- 
bonacci, and PAC algorithms. We show that practicable 
versions of mPAC and mFib (with small to) reach signif- 
icant cooling levels for 5-11 spin systems; Moreover, in 
the case of a 5-spin system, semi-optimal cooling levels 
(half the optimal cooling) are attained with reasonable 
run-time. Section V provides some analysis of SOPAC in 
case one would like to purify qubits as much as needed 
in order to obtain scalable NMRQC. We examine the 
resources required to obtain a highly polarized spin. 



II. SEMI-OPTIMAL PAC-BASED COOLING 
ALGORITHMS 

The Fernandez and Fibonacci algorithms described 
above illustrate the improved level of cooling attainable 
(asymptotically) by repeated compressions that involve 
the target spin. Practicable cooling algorithms that per- 
form a small number of such cycles (at each recursive 
level) offer reasonable cooling at a reasonable run-time. 
In this section we describe mPAC, a new PAC-based algo- 
rithm, which approaches nearly half the optimal cooling 
level with only a small number of cycles. 

The mPAC algorithm that we now present is a gen- 
eralization of PAC2 [26] (see Appendix A for details on 
PAC2). j e {1 ... J} is the purification level, k — 2j + 1 
is the bit index, B(k) stands for 3B-Comp on spins 
(k, k — 1, k — 2) which increases the bias of k, PT(A — > B) 
denotes a polarization transfer from bit A to bit B (or for 
simplicity, just a SWAP between their states), and Mq(1) 
is a reset, setting the bias of the reset spin, k = 1, to eo. 
The algorithm only uses elementary gates: a single gate 
operates either on a single spin or on a pair of spins (e.g., 
PT between adjacent spins), or on three spins (3B-Comp 
and PT between next-to-nearest neighbors). Mj(k) takes 
k = 2j + 1 spins at equilibrium and attributes bit k with 
a bias of Sj . 

mPAC: 



Mj(k) = [B(k)Mj-i(k - 2) PT(k - 2 -> k - 1) 



(2) 



Mj_i(fc - 2) ] m PT(k - 2 jfe) Mj-x{k - 2). 



With three spins 

Afi(3) - [B(3)M„(1)PT(1 ■ 

With five spins 

M 2 (5) = [B(5)M 1 (3)PT(3 



2)A/o(l)] m PT(1^3)Afo(l) 



4)M!(3)] m PT(3 5)Mi(3). 



The recursive formulas above are written from right to 
left, such that the first step of Mj(k) is reset of the reset 
spin, Mq(1), followed by PT from spin 1 to spin 3, and 
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FIG. 1. Cooling factors (e/eo, where eo = 10 5 ) for small 
spin-systems after mPAC with various m (see Eq. 2). 



to repetitions of the 4-step sequence in square brackets 
that ends with 3B-Comp. 

Notice that the number of spins required by mPAC to 
achieve a purification level J is the same as for PAC 2: 
n = 2 J + 1; but ej now depends on the number of cy- 
cles, to: £j = (2 — 2~ m )%o. For to = 1 we get back 
the algorithm PAC2, where Sj = (3/2)% , thus 1PAC 
= PAC2; we often retain here the original name, PAC2. 
Asymptotically, for to — > oo we get £j = 2 j £q. For three 
spins, C,B, A and to > 1, this is the Fernandez algorithm 
described above, and for sufficiently large to, spin C will 
acquire a final bias of E\ — 2eq [15, 22]. This is calculated 
via E\ = (2e + £i)/2 => e\ = 2e - Similarly, for five spins 
and sufficiently large to, the final bias of the 5 th spin is 
e 2 = (2ei + e 2 )/2 e 2 = 2et = 4e . 

The choice of to has a strong influence on the run- 
time, but fortunately, the polarization enhancement also 
increases rapidly with to. For small spin-systems (up to 
about 10 spins), very small values of to (2-6) are suf- 
ficient. Figure 1 compares the cooling factors obtained 
by such mPAC variants up to 13 spins. Notably, 6PAC 
cools to a similar extent as oo-PAC. It is also evident 
that 2PAC cools significantly better than the single-cycle 
variant (1PAC = PAC2). 

The run-time of mPAC (neglecting PT and 3B-Comp 
steps) is 



T(j) = T(j - 1) + 2T(j - 1)to 
= (1 + 2m)T(j - 1) = (1 - 



(3) 



2m) 3 



in reset-time units [27]. If to 
the run-time [15] is T(j) = 3 J . 
teger to > 1 



= 1 (mPAC is PAC2) 
In general, for any in- 



e 3 =e 3 - 1 {2-2- m ) = eo(2-2- m y. 

For instance, for 21 spins, j = 10, to = 7, we get Sj/so — 
985 (as long as the final bias is still much smaller than 
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1), and for the same j, but to = 4, Sj/eo = 745. For 11 
spins, j = 5, and m = 4 yields Sj/eo = 27.3. 

In general, for 2 j Sq <C 1 and sufficiently large m, spin 
fc = 2j + 1 is attributed a bias of Ej = 2ej-i = 2%o. Wc 
denote the asymptotic case (m —¥ oo) by ooPAC. 

We now generalize mPAC by replacing the constant to 
by a vector to. This added flexibility, in which a different 
value of to is associated with each cooling level, was 
found to be beneficial in the analysis of the Fibonacci 
algorithm [21], as we explain in section III. We call this 
new algorithm mPAC; it is defined as follows: 
mPAC: 

Mj(k) = [B(fc)Mj_i(fc - 2)PT(k - 2 -> k - 1) (4) 
M,-_i(fc - 2)p' PT(fc - 2 -> jfe) M,-_i(fc - 2), 

where (as before) j is the purification level. Mo(fc) de- 
notes resetting the bias of spin fc to e - 
With three spins 

Mi (3) = [B(3)M (1)PT(1 -> 2)M (l)] mi PT(1 -> 3)M (1) 
With five spins 

M 2 (5) = [B(5)Mi(3)PT(3 -> 4)Mi(3)] m2 PT(3 ->• 5)Mi(3) 

In general, Mj (k) cools spin fc to a bias of Sj — 2£j_i(l — 
2 -TOj '), as Mj is equivalent to performing the Fernandez 
algorithm on bits k , k — 1 , fc — 2 , where bits fc — 1 and fc — 2 
have equal initial biases £j-i- Therefore (see details in 
Appendix B) 

e j =eoH(2-2-^). 

3 

It is important to mention that m,j could be chosen ac- 
cording to various criteria; for instance, it could depend 
on the total number of spins n (or on J), in addition to 
its dependence on j. The algorithm requires J\j ;(1 + 2mj) 
reset steps: 

T(j)=T(j-l) + 2T(j-l) mj (5) 

3 

= (l + 2m j )T(j-l) = J](l + 2m j ). 

i=l 

Hence, the order in which the Wj appear is irrelevant. 

III. PRACTICABLE VARIANTS OF THE 
FIBONACCI ALGORITHM 

The Fibonacci algorithm is exhaustive, in the sense 
that a very large number of 3B-Comp steps are performed 
at each recursive level. Practicable variants may be ob- 
tained by limiting the number of compressions, as we 
describe here. An algorithm that reaches a distance of 
5* from the Fibonacci series was defined as follows: 
5-Fibonacci[21]: Run J"(n, n) 

T(n, k) = [T(n, k - 1) B{k)] m ^ k T(n, k - 1), (6) 



where F(n, 2) is a RESET step on bits 1 and 2 (described 
by Mq(1)PT(1 — > 2)Mo(l) in section II), m n .k are chosen 
such that the bias of bit fc is at least Jfe(l — 5„)eo ; and 
Fu is the fc th Fibonacci number. Here we choose 5% = 
2 fc_n_1 , such that £™ = ^. This condition sets m n .k — 
n — k + 2 (see appendix C). 
For three spins Fibonacci is 

J-(3, 3) = [J-(3, 2) B(3)] ro3 - 3 J-(3, 2) 
JT(3, 2) - M (l) PT(1 -> 2) M (l). 

This attributes spin number three with a bias of at least 
^3 (l — 2 _1 ) e = £o- For four spins 

J-(4,4) = [J"(4,3) S(4)] m4 - 4 J-(4,3) 
J"(4,3) = [J"(4,2) S(3)] m ^J"(4,2) 
J-(4,2) = M (1) PT(1 2) M (l), 

which attributes spins three and four with biases of at 
least J" 3 (1 - 2- 2 ) e = f £o, and J" 4 (l - 2" 1 ) e = f £o, 
respectively. 

The term <5-Fib hereinafter always refers to 5% = 
2*!-n-i ) guc j 1 that m n) fc = n — fc + 2. For large spin 
systems, say n > 8 or so, 5-Fib is not practicable, as it 
requires many cycles in the lower recursion levels. To cir- 
cumvent this problem, we fix the number of compression 
steps, such that m n ,k = m, and denote this variant of 
Fibonacci by mFib: 

mFib: Run J"(n, n) 

F{n, fc) = [J"(n, fc - 1) B{k)] m F(n, fc - 1). (7) 

Where J"(n, 2) is a RESET step on bits 1 and 2 as before, 
and m is any integer m > 1. With three spins 

J-(3, 3) = [-F(3, 2) B(3)] m T(3,2), 

and spin three is attributed a bias of 

e 3 = (l-2- m )( ei +e 2 )= 2(1-2-™). 

With four spins 

J-(4,4) = [^(4,3) £(4)] m J-(4,3) 
J-(4,3) = [^(4,2) S(3)] m J"(4,2), 

£3 is the same as for three spins, e 3 = 2(1 — 2~ m ), and 
spin four is attributed a bias of: 

e 4 = (l-2- m )(e 2 +e 3 ) 

= (l + 2(l-2- m ))(l-2- m ). 

In general, for n spins, the bias of the MSB is given by 
the recursive formula: 

e n = (l-2- m )( £n _ 1 +e„_ 2 ), 

where S\ — e 2 — 1. Specifically, we focus on cases where 
to is a small constant (3 — 5). For to = 2, it can be shown 
that mFib is outperformed by PAC3 [28] . 
The run-time of mFib is 

T(n,m) = (m + l)T(n- 1) (8) 
= (to + 1)"- 2 T(2) = 2(m + l) n - 2 . 
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IV. COMPARISON BETWEEN COOLING 
ALGORITHMS 

It is interesting to compare SOPAC to other algo- 
rithms. We first consider the cooling levels attained by 
each algorithm. ooPAC attributes to spin n = 2 J + 1 a 
bias of 2 J Eo. In comparison, with n spins, the PPA and 
all-bonacci reach 2 n_2 £o [29], and Fibonacci asymptot- 
ically approaches T n eo, where T n is the n th element of 
the Fibonacci series. While the PPA and the Fibonacci 
algorithm cool the entire spin-system, mPAC is defined 
so that it only polarizes the most significant bit (MSB). 
For a fair comparison of run-time we need to cool the 
entire string in mPAC as well. To accomplish this, suc- 
cessive applications of mPAC cool the less significant bits. 
Namely, the process 

M {1)PT{1 -I- 2)Mo(l)Mi(3)PT(3 4)Mi(3), . . . Mj{n) 

yields the asymptotic biases {2 J , 2 J_1 , 2 J ~ 1 
,8,8,4,4,2,2,1,1}. Consider the application of cooling 
algorithms to cool all spins (with initial biases of 0), until 
the biases are sufficiently close to the asymptotic biases; 
the resulting biases of the first seven spins (as long as e/ 
of the coldest spin is still much smaller than 1) are given 
here for ooPAC and other exhaustive algorithms, as well 
as for practicable and SOPAC algorithms: 

• ooPAC {2 J ,..., 8, 4, 4, 2, 2, 1,1} 

• PPA and all-bonacci {2"- 2 , . . . , 32, 16, 8, 4, 2, 1, 1} 

• Fibonacci {J"„, . . . , 13, 8, 5, 3, 2, 1, 1} 

. PAC2= 1PAC {(3/2) J ...,3|,2i,2i,li, 11,1,1} 
. PAC3 {...,5±,3§,2§, 2,11,1,1} 

• 2PAC {. . . , 5.36, 3.06, 3.06, 1.75, 1.75, 1, 1} 

• 4PAC {. . . , 7.27, 3.75, 3.75, 1.94, 1.94, 1, 1} 

• 6PAC {. . . , 7.81, 3.94, 3.94, 1.98, 1.98, 1, 1} 

• 3Fib {. . . , 7.81, 5.29, 3.64, 2.41, 1.75, 1, 1} 

• 4Fib {. . . , 10.2, 6.54, 4.28, 2.70, 1.88, 1, 1} 

The bias configurations are given in units of the initial 
bias, £q. 5-Fib is not included, as for each total number 
of spins it produces different bias series. For n = 7, <5-Fib 
yields the biases {8.27, 6.49, 4.54, 2.88, 1.97, 1, 1}. 

We consider a small spin-system comprising five spins. 
Table I compares the biases (for the MSB) obtained by 
previous algorithms (top), mPAC (middle section), and 
mFib (bottom), as well as the number of resets required 
to create the entire bias series [30]. Note that the PPA 
cools better when more resets are allowed, approaching 
the limit of 8£o for AC with 5 spins. With only 28 resets, 
the PPA attains a semi-optimal cooling level of 4e ; For 
a similar cooling level, 4PAC and 6PAC require 101 and 
197 resets, respectively. After 99 resets, the PPA obtains 



a near-optimal cooling level of 7eq. Table II presents a 
similar comparison for 7 spins. Table III and Table IV 
present similar comparisons for 9 and 11 spins, respec- 
tively; a spin-system of comparable size was recently used 
for benchmarking quantum control methods [4]. Ap- 
pendix D compares the number of spins and run-time 
required by each algorithm to achieve several small cool- 
ing factors. 

The run-time analysis of the PAC algorithms is conve- 
niently expressed as a function of the purification level, 
J (see Appendix A). The entire spin-system is cooled 
by 2 J — 1 successive applications of each algorithm, as 
shown above for mPAC. In contrast, the Fibonacci algo- 
rithm and the PPA were designed to generate the entire 
series of biases. The run-time of <5-Fib is given by n! (see 
calculation in appendix C), and the run-time of the PPA 
was obtained by a computer simulation that iterates be- 
tween the two steps of the algorithm. 





bias (eo) 


Run-time 
(reset steps) 


5-Fib 


3.29 


120 


PAC2 


2.25 


17 


PAC3 


2.75 


29 


PPA 


4.03 


28 


PPA 


7.00 


99 


2PAC 


3.06 


37 


4PAC 


3.75 


101 


6PAC 


3.94 


197 


3Fib 


3.64 


128 


4Fib 


4.28 


250 



TABLE I. Cooling a 5-spin system by various algorithms. The 
optimal cooling level is 8eo. 





bias (eo) 


Run-time 
(reset steps) 


5-Fib 


8.27 


5040 


PAC2 


3.38 


53 


PAC3 


5.13 


169 


PPA 


8.02 


104 


PPA 


16.0 


428 


2PAC 


5.36 


187 


4PAC 


7.27 


911 


6PAC 


7.81 


2563 


3Fib 


7.81 


2048 


4Fib 


10.15 


6250 



TABLE II. Cooling a 7-spin system by various algorithms. 
The optimal cooling level is 32eo. 
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KlQO f C j-* 1 

UlcLo ^c() J 


Run-time 
(reset steps) 


(5-Fib 


21.5 


362880 


PAC2 


5.06 


161 


PAC3 


9.56 


985 


PPA 


32.0 


1639 


PPA 


64.0 


6836 


2PAC 


9.38 


937 


4PAC 


14.1 


8201 


6PAC 


15.5 


33321 


3Fib 


16.9 


32768 


4Fib 


24.2 


156250 



eration: 



TABLE III. Cooling a 9-spin system by various algorithms. 
The optimal cooling level is 128eo- 





bias (eo) 


Run-time 
(reset steps) 


5-Fib 


56.0 


39916800 


PAC2 


7.59 


485 


PAC3 


17.8 


5741 


PPA 


64.0 


6456 


PPA 


256 


109323 


2PAC 


16.4 


4687 


4PAC 


27.3 


73811 


6PAC 


30.8 


433175 


3Fib 


36.4 


524288 


4Fib 


57.7 


3906250 



TABLE IV. Cooling an 11-spin system by various algorithms. 
The optimal cooling level is 512so. 



V. PURIFICATION OF QUBITS 

The practicable nature of the algorithmic cooling de- 
scribed in the previous sections might hold significant po- 
tential for deriving scalable NMR quantum computation. 
The pseudopure states used in NMRQC were found to 
scale well for sufficiently pure spins, for which e = 1 — 2(5, 
where the error probability, 6, satisfies 5 <C 1, see [8]. For 
small spin systems, <5 <~ 0.2 is already useful [31], while 
larger systems require (6 of the order of 1/n) [8] . Evalu- 
ation of cooling algorithms for qubit purification requires 
exact calculations that go beyond the approximation of 
e <C 1, which was used so far. Here we describe the exact 
result of 3B-Comp, and use it in the general analysis of 
mPAC. 

The precise result of 3B-Comp is given by [15, 22]: 

/ £3 + e 2 + £l — £3 £ 2£l 



Performing 3B-Comp iteratively while replenishing e\ 
and £2 after each compression will give, in the m th it- 



Jrn) 



_(m- 



>-i), „ 

■■l-£ 3 £2£l 



£2+£i+4" 1 1) (l~£2£i) 



In the case of mPAC, ei = £2 and therefore 



Defining 



_(m) _ £3 I 1 £ l) 



1 - F 2 

A(e) = 



+ £1. 



2 ' 



(9) 



(10) 



we solve Eq. 9. 



A{ £l ) 2 s { ™- 2) + {A{s l ) + l)s l 

m—1 



i=0 



= A{e l Tef ) +e 1 Y J A{e 1 ) 
In our scheme, — e\ and so 

m 

ei^^(ei) 



1+m 



£l 



1-A( £l ) 
l-A( £l ) 



(11) 



which can be generalized to the following recursive for- 
mula: 



_(m) 



e k -2 22A(e k - 2 ) 1 =£ 

i=0 



(m) 
fe-2" 



-A(etl) 



l+m 



I- A 



(4- } 2 ) ' 



(12) 

Cooling curves for several mPAC variants are pre- 
sented in Figure 2, starting from a typical electron polar- 
ization of 1%. The curve for 6PAC is similar to oo-PAC 
(see Figure 4 in Appendix E for a close-up view); 6PAC 
requires about 13 spins to attain a semi-optimal bias, and 
with about 17 spins it reaches polarization near unity. 
Starting from hyperpolarization (eo = 10%), 6PAC ex- 
ceeds 50% polarization with only 7 spins (see Figure 3 
and Figure 5); In this case, 11 spins are sufficient to pu- 
rify one spin to a high degree. 

Table V gives the resources required by SOPAC in 
order to reach the threshold of scalable NMRQC, e = 
0.6 [8]; The left pane assumes an initial polarization of 
1%, a typical electron polarization, while the right pane 
assumes eo = 10%. Results for mFib were obtained nu- 
merically. Note that for m — 6, mPAC requires the same 
amount of spins and more reset steps than 4PAC, and 
produces slightly higher biases. 



G 



1 
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FIG. 2. Polarization biases (e2j+i) for mPAC with various m 
as a function of purification level j (defined in Eq. 2). The 
initial bias is £q = 0.01. The curve for ooPAC (i.e., m — ¥ oo) 
was calculated from Eq. B2. 




2 3 4 5 

purification level (j) 

FIG. 3. Polarization biases (£2j+i) for practicable mPAC vari- 
ants as a function of purification level j. The initial bias is 

£0=0.1. 



Fault-tolerant quantum computing 

For the purpose of fault tolerant quantum computing, 
it is essential to obtain highly pure spins, where the error 
probability, S — -jp, is well below 1%. Table VI shows 
the number of resets and spins required by mPAC vari- 
ants to achieve e > 0.9999 for the MSB, starting from 
£o = 0.1, which is a typical electron spin polarization 
at low temperature. For m — 6, mPAC still requires 13 
spins and more reset steps than 4PAC (the bias achieved 
by 6PAC is 0.99999). 

It is notable that 2PAC and 4PAC achieve such high 





£ > 0.6 


£ = 0.01 


£0 =0.1 


alg. 


#spins 


run-time 


bias 


#spins 


run-time 


bias 


1PAC 


23 


354293 


0.72 


11 


485 


0.66 


2PAC 


17 


585937 


0.72 


9 


937 


0.76 


4PAC 


15 


5978711 


0.78 


7 


911 


0.63 


3Fib 


13 


97656250 


0.72 


7 


2048 


0.67 


4Fib 


12 


19531250 


0.72 


7 


6250 


0.78 



TABLE V. Number of spins and run-time required by SOPAC 
to purify the MSB to a polarization of at least 60%, starting 
from an initial polarization of 1% (left) or 10% (right). For the 
purification from 1%, 3PAC requires 15 spins, like 4PAC, and 
less resets (1098057), but reaches a lower purification (68%). 



purity with only six more spins (three additional purifi- 
cation levels) with respect to the minimal purification 
(~ 60% - 80% purity) obtained in Table V. This rapid 
purification stems from the behavior of mPAC at low er- 
ror probabilities, which may be obtained by substituting 
e by 1 — 25 in Eq. 12, and considering the first three terms 
of the summation (i = 0, 1, 2): 

1 - 25 j+l = (1 - 25j) [l + 25,(1 - 5j) + 4^(1 - S 3 ) 2 ] , 

where Sj is the error probability at purification level j, 
and A(e) = A(l - 25) = 26(1 - 5) (see Eq. 10). At 



low error probabilities, 5j+i 



5? As a result of 



this rapid convergence, the number of spins required by 
mPAC to reach 99.9% purity is sufficient to produce a 
markedly higher purity of 99.99% (or even 99.999% for 
some variants). 

Table VI includes two mFib SOPAC variants, which 
require fewer spins than the practicable mPAC variants; 
Both 3Fib and 4Fib require 11 spins, while mPAC re- 
quires at least 13 spins. The run-times and cooling levels 
of 3Fib and 4PAC are similar, while 4Fib has a much 
higher run-time and reaches extremely high purification 
(99.9999%). 

Obviously, only the case of purifying from e = 0.1 to 
above 60% polarization can be considered as a potential 
near future application. The other cases presented here 
require an unreasonable run-time. 



VI. DISCUSSION 

We have introduced SOPAC, semi-optimal practicable 
cooling algorithms that achieve better cooling than the 
PAC algorithms by performing a small number of cy- 
cles at each recursive level. The mPAC and its special 
instance mPAC (constant m) offer improved cooling for 
small spin systems at a reasonable increase in run-time. 
Similar results were obtained by means of the mFib al- 
gorithm, which requires a constant number (m) of 3-bit 
compressions at each recursive level. We compared the 
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alg. 


#spins 


run-time 


final bias 


1PAC 


19 


39365 


0.999996 


2PAC 


15 


117187 


0.999999 


4PAC 


13 


664301 


0.999984 


3Fib 


11 


524288 


0.999877 


4Fib 


11 


3906250 


0.999999 



TABLE VI. Number of spins and run-time required by 
SOPAC to achieve a minimum bias of 0.9999, starting from 
£o = 0.1. With the same number of spins as 4PAC, 3PAC 
reaches a slightly lower bias (0.999938) using fewer resets 
(156865). 



performance of mPAC and mFib (with small m) with 
the previous algorithms - the PAC algorithms (PAC2 and 
PAC3), the PPA, and <S-Fib. While the PPA is proven 
to be optimal, it is not expected to be the first choice in 
practice, since an efficient implementation of its SORT 
step might not be found. The SOPAC algorithms allow 
flexible tuning of the degree of cooling, thereby bridging 
the large gap between the elementary cooling of PAC, 
and the extensive, yet unattainable, cooling of exhaus- 
tive algorithms. 

AC was mentioned as an enabling technique for various 
quantum computing schemes [32-35]. The simplicity of 
implementation and semi-optimal nature of SOPAC ren- 
ders it potentially useful for deriving spin-based quantum 
computers in the more distant future [8, 32-35]. To this 
end, the new practicable cooling algorithms presented 
here may yield sufficiently pure qubits, starting from rea- 
sonably small spin systems initially polarized by nearby 
electrons at room temperature (e = 0.01) or at low tem- 
perature (eo = 0.1). In the second case, both 3Fib and 
4PAC require only seven spins to cool the MSB beyond 
60% polarization, permitting scalable NMRQC [8]. With 
about five additional spins, these algorithms can ideally 
produce highly purified spins (e > 0.9999) that are suit- 
able for fault tolerant quantum computing. While a sim- 
ilar 12-spin system was recently manipulated in liquid- 
state NMR, more spins would be required in practice due 
to various factors, such as actual thermalization times 
and decoherence in between reset steps due to spin-spin 
relaxation. 

AC may contribute to NMR spectral editing, wherein 
part of the spectrum that corresponds to a particular spin 



of interest is selectively enhanced with respect to overlap- 
ping signals. In such cases, SOPAC can potentially offer 
selective enhancement much beyond simple polarization 
transfer. 

Last but not least, in real life, the reset spins do not re- 
lax infinitely faster than the computation spins (see [36- 
41] for viewing AC as a novel type of heat-engine). Let TZ 
denote the ratio between the relaxation time of the com- 
putation spins and the relaxation time of the reset spins, 
TZ = Ti(comp.)/Ti (reset). Ref [28] thoroughly analyzes 
various algorithms, including 2PAC suggested here, for 
several values of TZ, and found (as expected) that when 
the ratio is much larger than the run-time (number of 
reset steps) required by the algorithm, TZ T run _ t i me , 
the results here still apply with only minor corrections. 
However, when 1Z is in the range of about T rtm _ time to 
10T run -time, the deviations become significant. For ex- 
ample, Ref [28] considers the application of 2PAC to 7 
spins, starting from £o <C 1. When 1Z — 10000 and the 
duration of each reset step is 5Ti (reset), cooling all the 
spins attributes the MSB with a final bias of 5.11 (in units 
of £o) ; slightly below the ideal bias (5.36, see Table II); 
lower ratios, 1Z — 1000 and 1Z = 100, significantly re- 
duce the final bias, to 3.63 and 1.07, respectively. When 
only the MSB is cooled, the same bias, 5.36, is reached 
in the ideal case, while the run-time is shortened from 
187 to 125 reset steps, hence significantly higher biases 
are obtained; for 1Z = 10000 (TZ = 1000, TZ = 100) the 
resulting biases arc 5.24 (respectively 4.58,2.56). The 
final bias may be further increased by optimizing the du- 
ration of the reset step (in particular for small TZ); For 
example, when TZ = 100, reducing the reset duration to 
I.8T1 (reset) increases the final bias from 2.56 to 3.05, 
when cooling only the MSB, while reducing the reset pe- 
riod further, to 0.91Ti (reset), nearly doubles the final 
bias, from 1.07 to 1.97, when all spins are cooled (see [28] 
for more details). Notably, TZ <~ 10000 was achieved in 
solid-state NMR, where spin-diffusion was used for rapid 
repolarization [14, 16]. 
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Appendices 



Appendix B: Analysis of mPAC 



Appendix A: Formal descriptions of PAC algorithms 

PAC2 on 2 J + 1 spins can reach a purification level 
J on the MSB spin (k = 2 J + 1) [9]. The procedure 
M (l) is defined as a reset step on spin k = 1. The 
procedure Mj(k) is defined recursively to cool the k th 
spin to a purification level j € 1 ... J (i.e. £j = (3/2) J ). 
i3(&:) denotes 3B-Comp on spins (k, k — l,k — 2) which 
increases the bias of k from level j — 1 to level j, and 
PT(A — > P) denotes a polarization transfer from bit A 
to bit P. 

Practicable algorithmic cooling 2 (PAC2): 
For j €{1,..., J} 

Mj-(jfe) = B(fc)M,-_i(ft - 2)PT(k - 2 -> jfc - 1) 

Mj-_i(fe - 2)PT(fc- 2 -> fe)Mj-_i(A; - 2), 

With Mo(l) being reset on spin 1. For example, to cool 
spin 3 to the first level of purification we use: 

Mi(3)=B(3)M (l) 

PT(1 -> 2)M (1)PT(1 3)M (1), 

and to cool spin 5 to the second purification level we use 

M 2 (5) =B(5)Mi(3) 

PT(3 -> 4)Mi(3)PT(3 5)Mi(3). 

PAC2 requires J = [logj 5 (£j/£o)] , a s every 3B-Comp 
step boosts the bias (ej <C 1) by a factor of 1.5 [9]. 
The calculated run-time of PAC2 is 3 J , for cooling the 
MSB [9]. 

PAC3 is described as follows [15]: 
Practicable algorithmic cooling 3 (PAC3): 
for fee {2, . . . , J} (the reset spin is at index 1) 

M k (k+1) =B(& + l)M fc _ 2 (fc-l) 

M k ^(k)PT(k -> k + l)M fc _i(A!), 

with M (l) = RESET (of the reset spin), and 

Mi (2) = PT(1 — > 2)M (1). For three spins PAC3 is 

expressed as 

M 2 (3) = S(3)M (1) Mi (2) PT(2 -> 3)Mi(2), 
and with four spins PAC3 is described by 

M 3 (4) = S(4)Mi(2) M 2 (3) PT(3 -> 4)M 2 (3). 
The resulting bias series on 7 spins is 



5i 3^ 2^ 2 2 1 ] * 

8' 4' 4' '2' ' 



The run-time of PAC3 is given by the recursive for- 
mula [15]: Tj = 2T 7 _ 1 + Tj_ 2 , where T = T 1 = 1, 
and J = n — 1. 



We can generalize the expression for the achieved bias 
(Eq. 12) to 



'fc-2 



1 



^fc-2 J 



1+mj 



^fc-2 ) 



The limit of this algorithm is given by 



lim 



-(oo) 
-fe-2 



.4 



(4-0 



(Bl) 



(B2) 



.(oo) 



where ef"' = = Co- 
in the limit of e k <C 1 we can use the first order 
approximation in e k . The bias assigned to the target 
bit, k, thus depends only on {m^i = 1, . . . j} (i.e., with 
no significance to the order of appearance). For exam- 
ple, mPAC on five spins which uses m x = 3, m 2 = 5 
will set the coldest bit to the same bias as the sequence 
mi = 5, m 2 = 3. To prove this let us first notice that 



Eq. 10 becomes A(e) 
analytical solution 



£<1 



Eq. Bl then has a simple 



S m j) 



2*rr J ( i -2- 



(B3) 



i=l 



Appendix C: Derivation of the run-time of 5-Fib 

The choice of m n ,k = n — k + 2 in Eq. 6 yields biases 
that approach the Fibonacci series (J"„, T n -\, ■ ■ ■), up to 
a factor of 1 — S n ^ = \ — 2 k ~ n ~ 1 . We assume by induction 
that bits k — 1 and k — 2 attained the biases £o a i > 
£ Jfe-i(l - 2 k - n - 2 ) and e a 2 > £ Pfe- 2 (l - 2 fe -"- 3 ), 
respectively. We shall now prove that bit k will attain a 
bias e x > £o-Pfc(l - 2 fe ~"~ 1 ). 

The effect of the 3B-Comp gate applied to the bits 
k — 2, k — 1, k with the corresponding biases £o»i, £o a 2, 
and eox, assigns bit k the bias Eqx — > £o(ai + a 2 + x)/2 
(with ai£o, a 2 £o <C 1). In each iteration within the fc th re- 
cursive level, x converges towards the unique fixed point 
of this transformation, x = a\ + a 2 . The convergence is 
rapid; after m n _ k iterations, starting at x = 0, Fibonacci 
achieves x — (a\ + a 2 )(l - 2~ m ™. fe ). Note that 

oi + a 2 > J- fe _i(l - 2 fe -"- 2 ) + 7fc_ 2 (l - 2 fc -"- 3 ) 
= T k — 2 k n 3 (2J 7 k-i + Tk--i) 



T k 



-ik—n—3 



Pfe+i > -Pfc(l — 2 



fc-ri-2 



), 



where the last inequality is due to .Ffe+i < 27^. There- 
fore, by choosing m n . k = n — k + 2, the Fibonacci algo- 
rithm will achieve 



x> T k (l-2 k - n - z )(l 
>T k (l-2 k - n - 1 ), 



')=T k (l-2 



k-n-1\1 
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where the last line is due to (1 — a) 2 > 1 — 2a, for any 
positive a. □ 



Appendix D: Comparison of algorithms for small 
target cooling factors 



According to the definition of Fibonacci in Eq. 6, the 
run-time T(n, k) of the operation J"(n, k) is given by 



T(n,k) = {m nik + l)T(n,k-l) 

k k 

= T(n, 2) J] (m„, fc + 1) = 2 JJ (n - j + 3), 

i=3 j=3 



where only the run-time of RESET steps is taken into 
account. T(n, 2) = 2 because J"(n, 2) involves two RE- 
SET steps. Therefore, the run-time of the full Fibonacci 
algorithm, ^(n, n), is: 



n n 

T(n,n) = 2 JJ(n- j + 3) =2 J[i = n\. 

j=3 i=3 



Only half the run-time is required if the two reset steps 
in J"(n, 2) are performed in parallel. 



It is interesting to compare the performance of SOPAC 
with PAC and other cooling algorithms for fixed target 
cooling factors. Tables VII, VIII, IX, and X list the total 
number of spins and reset steps required by each cooling 
algorithm to reach biases of at least 3eo, 7eo, lleo, and 
15£o> respectively. 





#spins 


Run-time 
(reset steps) 


5-Fib 


5 


120 


PAC 2 


7 


53 


PAC3 


6 


70 


PPA 


4 


16 


2PAC 


5 


37 


3Fib 


5 


128 



TABLE VII. A comparison of various algorithms for a target 
bias of at least 3eo- The only SOPAC variants shown are 
2PAC and 3Fib, as other practicable variants require the same 
amount of spins and more reset steps. 



Appendix E: Close-up view of cooling curves for 
4PAC, 6PAC and ooPAC 
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Tropins 


Run-time 
(reset steps) 


5-Fib 


7 


5040 


PAC2 


11 


485 


PAC3 


8 


408 


PPA 


5 


97 


2PAC 


9 


937 


4PAC 


7 


911 


6PAC 


7 


2563 


3Fib 


7 


2048 


4Fib 


7 


6250 



TABLE VIII. Performance of various algorithms for a target 
bias of at least 7eo • 





#spins 


Run-time 
(reset steps) 


5-Fib 


8 


40320 


PAC2 


13 


1457 


PAC3 


10 


2378 


PPA 


6 


204 


2PAC 


11 


4687 


4PAC 


9 


8201 


6PAC 


9 


33321 


3Fib 


8 


8192 


4Fib 


8 


31250 



TABLE IX. Performance of various algorithms for a target 
bias of at least lleo- 





#spins 


Run-time 
(reset steps) 


5-Fib 


9 


362880 


PAC2 


15 


4373 


PAC3 


11 


5741 


PPA 


6 


529 


2PAC 


11 


4687 


4PAC 


11 


73811 


6PAC 


9 


33321 


3Fib 


9 


32768 


4Fib 


8 


31250 



TABLE X. Performance of different algorithms for a target 
bias of at least 15eo- 
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FIG. 4. Close-up view of the polarization biases (E2j+i) for 
4PAC, 6PAC and ooPAC as a function of the purification level 
j, starting from an initial bias of so = 0.01. For more details 
see Figure 2. 
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FIG. 5. Close-up view of the polarization biases (e2j+i) for 
4PAC, 6PAC and ooPAC as a function of the purification level 
j, starting from an initial bias of eq = 0.1. For more details 
see Figure 3. 
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